Artifact of the phonon-induced localization by variational calculations in the 

spin-boson model 
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We present energy and free energy analyses on all variational schemes used in the spin-boson 
model at both T = and T 7^ 0. It is found that all the variational schemes have fail points, at 
where the variational schemes fail to provide a lower energy (or a lower free energy at T 7^ 0) than the 
displaced-oscillator ground state and therefore the variational ground state becomes unstable, which 
results in a transition from a variational ground state to a displaced oscillator ground state when 
the fail point is reached. Such transitions are always misidentified as crossover from a delocalized 
to localized phases in variational calculations, leading to an artifact of phonon-induced localization. 
Physics origin of the fail points and explanations for different transition behaviors with different 
spectral functions are found by studying the fail points of the variational schemes in the single mode 
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I. INTRODUCTION 



Variational method is one of the important tools in 
tracing the ground state of many body systemSi^ In deal- 
ing with systems having interaction with phonons, the 
variational method based on displaced-oscillator state 
(DOS) plays an important role. Some well-known 
examples in condensed matter physics are polaron- 
phonon, exciton-phonon systems j-^j^s^ and spin-boson 
model, which is an important toy model in dissipative 
quantum systems In fact, the coupling to a phonon 
bath in all the mentioned systems has the same bilin- 
ear coupling form. One important issue in the men- 
tioned systems is the so-called phonon-induced local- 
ization. In exciton-phonon (or polaron-phonon) sys- 
tems, the phonon-induced localization is called the self- 
trapped transition, while in the spin-boson model, it is 
stated as the crossover from a delocalized to localized 
phasesi^'i'^ Such a crossover is now considered as some 
kind of quantum phase transition called boundary phase 
transition^iiS In general, variational calculation has been 
serving as an important tool in dealing with the phonon- 
induced localization, however, the predicted transitions 
in exciton-phonon (or polaron-phonon) systems turned 
out to be false ones but the answer to why the variational 
calculation fails was not found. Variational calculation 
on the crossover problem in the spin-boson model has 
a long history and the first variational scheme is based 
on the DOS.-i^ii^'i^It was shown that, in strong cou- 
pling regime, the coupling to phonon bath can lead to 
both displacement and deformation of the oscillator and 
the so-called displaced-squeezed-state (DSS) is energy fa- 
vorite, leading to a variational calculation based on the 
DSS.^^'^^ Later, a hybrid variational scheme based on the 
DSS was also developedJ^ At finite temperature regime, 
the variational calculation based on the DOS was done 
for both the Ohmic and sub-Ohmic cases and a discon- 
tinuous crossover was found )^'^i^^'^^i^° and recently, it is 
shown that the variational calculation based on the DSS 



fails in thcrmodynamical limit. 

Despite of its long history, the result from variational 
calculations is far from satisfactory. The basic prob- 
lem is that the predicted crossover point is variational- 
scheme dependent. For example, in the Ohmic case, the 
crossover point by the variational calculation based on 
the DOS was ac = Ir^ while that based on the DSS gave 
ttc = 2j^ii^ and the hybrid variational scheme predicted 
a ac lay between. 1^ The transition behavior predicted 
by the variational calculations is also in controversy. 
The variational calculations based on both the DOS and 
DSS predicted a continuous transition in the Ohmic case 
but discontinuous transition behavior was also found in 
the hybrid variational schemeJ^ In the sub-Ohmic case, 
all the variational calculations predicted a discontinuous 
transitioufiSii^i but non-perturbative numerical calcula- 
tions, like numerical renormalization group (NRG) and 
density- matrix renormalization group (DMRG), found a 
continuous transition insteadi^i^iiS^ The situation in the 
case of T 7^ is also confusing. Variational calculations 
predicted a discontinuous transition at T 7^ even in 
the super-Ohmic case?^ a result which is in confliction 
with the known conclusion that the crossover only hap- 
pens in the Ohmic and sub-Ohmic casesj^siiS All these 
discrepancies make the crossover predicted by the vari- 
ational calculations in question, e.g., is it an artifact as 
in the case of polaron-phonon systems? Clearing up this 
issue is important for understanding the result by varia- 
tional calculations and further application of variational 
methods to other systems. This is the motivation of the 
present paper. 

By employing energy analysis on the variational 
schemes used in the spin-boson model, we shown that 
the predicted crossover is just an artifact of the varia- 
tional schemes. We will show that all the variational 
schemes mentioned above have fail points, at where vari- 
ational schemes fail to provide a lower energy than that 
of the displaced-oscillator ground state, leading to a tran- 
sition from the variational ground state to the displaced- 
oscillator ground state when the fail point is reached. 
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In variational calculations, such a transition is always 
misidentified as a crossover from a delocalized to localized 
phases and the fail point is misinterpreted as a crossover 
point. The result we found can help to resolve the con- 
troversy mentioned above. 

The rest of the paper is organized as follows: The fail 
points of the variational schemes at both T = and 
T 7^ are demonstrated in the next section. In Sec III, 
physical origin of the fail point and different transition 
behaviors in both the Ohmic and sub-Ohmic cases are 
analyzed by studying the variational schemes in the single 
mode case. Conclusion and discussion are presented in 
the last section. 



II. FAIL POINTS OF THE VARIATIONAL 
SCHEMES USED IN THE SPIN-BOSON MODEL 

The Hamiltonian of the spin-boson model is given by 
(setting h=l)^ 

k k 

where ai (i = x,y,z) is the Pauli matrix, bk (bl) is the 
annihilation (creation) operator of the fcth phonon mode 
with energy ivk and Ck is the corresponding coupling pa- 
rameter. As is known, the solution of this model depends 
on the so-called bath spectral function that is defined as 
J(ti>) ~ ^'l2k^k^i^ ~ Usually the spectral func- 

tion has a power law form, i.e.. 



J(w) = -acj'tj^ 



< uj < ujd, 



(2) 



where uod is the cut-off frequency and a is a dimension- 
less coupling strength which characterizes the dissipation 
strength. ^'^ The property of the bath is characterized by 
the parameter s, i.e., 0<s<l,s = l, and s > 1 
represent respectively the sub-Ohmic, Ohmic, and super- 
Ohmic cases. 

It is well known that, with the bilinear coupling form 
given in the last term of Eq.(l), the Hamiltonian can be 
diagonalized by the so-called displaced-oscillator trans- 
formation in the case of Aq = 0. The ground state in 
this case is the displaced-oscillator ground state^i^i^^ 



l<fo> 



1 



x/2 



(l</'+>IT> + l0->U», 



(3) 



where |(/)±) = exp{=F I]fc(cfe/wfc) (6fc-6^)}|0), is the DOS, 
while I t) and | [) are the eigenstates of CTz. All the 
variational schemes used in the spin-boson model start 
from the DOS; analysis in details will be presented in 
the following. 



A. Fail points of the variational schemes at T = 

As we have mentioned in the Sec I, there are three 
kinds of variational schemes used in the spin-boson model 



at r = and the key difference is the trial ground 
state. The trial ground state for the variational calcu- 
lation based on the DOS is chosen as^^ 



1 



where 



|*o)-^(l^0-f)IT> + IV^-)U)), 



l^±) =exp|T^5fe (^fe-^I.)| |0), 



(4) 



(5) 



and gk is the variational parameter. The trial ground 
state based on the DSS is given byi^ii^ 



where 



|*o>-^(l¥'+>IT> + l^-)U)), 



(6) 



l^±> =exp{T^- (fo.-5i.)}exp{-^7fe(fo?.-&r)}|0), 

k ^ k 

(7) 

is the DSS and 7^ is the variational parameter. The third 
variational scheme is a hybrid scheme of the above two 
variational schemes and the trial ground state is^ 

1 



with 



|$o) = ^(|^+)|T) + l'^-)U)), 



exp{T X! 9k{bk -bl)} cxp{- ^ jk{ 

k k 



(8) 
(9) 

here both gk and jk are the variational parameters. 

All the variational schemes have almost the same calcu- 
lation process and the main routine is as follows. Firstly, 
one used the trial ground state \^tgs) to calculate the 
trial ground state energy Et - 
following variational condition 



<Ptgs\H\^tgs), then the 



SEt/Sgk - SEt/Sjk = 0, 



(10) 



can help to determine the variational parameters. In 
what follows, the trial ground state with the variational 
parameters determined from the above variational condi- 
tion is called the variational ground state {"^v)- Variation 
of the tunneling splitting of the variational ground state 
with the dissipative strength a can tell how the crossover 
happens. Since all three variational calculations give the 
same crossover behavior in both the Ohmic and sub- 
Ohmic cases, here we shall use the variational calculation 
based on the DOS as an example to show the details. 

For variational scheme based on the DOS, the varia- 
tional condition (10) leads to 



9k 



Ck 



ujk + Ao cxp{-2X;fc.9fc}' 



(11) 



from which the self-consistent equation of the dressing 
factor Ki can be foun d^^'^'^i^^ 



i^i =/(i^i) = cxp|-^^ 



KiAa)' 



;dLu\, (12) 
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here Ki = A/ Aq and A is the tunnehng spHtting in pres- 
ence of dissipation. The trivial and non-zero solutions 
of Eq. (12) represent respectively the localized (A = 0) 
and delocalized (A ^ 0) phases. Variational calculation 
determines the crossover by examining the evolution of 
solutions of Eq. (12) with the dissipation strength a. It 
turned out that, in the Ohmic case, there is only one 
non-zero solution which decreases to zero continuously 
as a approaches to 1 in the limit of Ao/cij£) <C 1, showing 
a continuous crossover happens at aj, = 1- The situa- 
tion in the sub-Ohmic case is more complicated; there 
are two non-zero solutions {K[ and K'() which disappear 
abruptly when a reaches some critical value. It has been 
shown that the crossover in this case is discontinuous 
and the crossover point should be determined according 
to the Landau theory.^^'^^ The situation for other vari- 
ational schemes is the same except for the values of the 
critical points. It should be noted that the discontinu- 
ous transition behavior in the Ohmic case found in the 
hybrid variational schemeii can also be found in other 
variational schemes and the key point is to increase the 
value of Ao/w_D. We find that, for Aq/loo > 0.5, a 
discontinuous transition appears in the Ohmic case for 
the variational scheme based on the DOS. The change of 
transition behaviors is shown in Fig. 1. Explanation for 
this result will be presented in next section. 
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FIG. 1: (Color online) The evolution of solutions of the self- 
consistent equation derived from the variational calculation 
based on the DOS in the Ohmic case. In the case of Aq/ujd ~ 
0.1, there is only one non-zero solution and the transition 
behavior is continuous. However, in the case of Aq/lod = 1, 
there are two non-zero solutions K[ and A"" and the transition 
behavior becomes discontinuous as a approaches ac- 

Now we turn to show how variational schemes fail as 
a increases. As we have mentioned, the starting point 
of all the variational schemes is DOS. Consequently, for 
a successful variational scheme, the variational ground 
states I'^v) should be a better approximation to the true 
ground state than the displaced-oscillator ground state 
|$o) 1 which implies that j^E't,) should have a lower energy 



than \^o), e.g., the following condition 



(13) 



should hold in a variational calculation, otherwise the 
variational scheme fails. Surprisingly, it is found that 
all the variational schemes fail when a increases to some 
critical value in both the Ohmic and sub-Ohmic cases. 
In the sub-Ohmic case, the variational ground state has 
a lower energy than the displaced-oscillator ground state 
when a <^ ac- As a increases, both Ey and Eq decreases 
but Eq decreases in a way faster than Ey , we have Ey = 
Eq when a = ac, and Ey > Eq as a > ac, showing that 
the variational scheme fails at a — ac- Typical results 
for the variational calculation based on the DOS in the 
case of s = 0.5 and Aq/ojd — 0.1 are shown in Fig. 
2. For the Ohmic case with Ao/ojd ^ 1, the energy 
difference Ey — Eq decreases to zero continuously as a — > 
ac- The results of the variational calculation based on the 
DOS is shown in Fig. 3. However, as Aq/wd increases to 
some larger values, the self-consistent equation has two 
non-zero solutions as shown in Fig. 1, and the situation 
becomes the same as in the sub-Ohmic case. 

Analytically, the fail point ac of the variational scheme 
can be determined by 



Ey — En. 



(14) 



In the case of the variational scheme based on the DOS, 
this leads to 



ac 



Ao(Xi-Ko) 



Jo ^ujd' \ui+KiAqJ 



(15) 



where Kq = (0+|(/)_) = exp{ — ;| J^" ^^-^du)}, is the adi- 
abatic dressing factor. In the Ohmic case, we have s = 1 
and Kq — and the fail point is 



ac ^1 + KiiAQ/LOo). 



(16) 



The above equation is just the crossover point predicted 
by the variational calculation in the limit of Aq /lod <^ 
Ifl^ili On the other hand, for the variational scheme 
based on the DSS, the fail point equation tells 



^AqD = ^sinh^(27fc)wfe| 



(17) 



where is determined by e'^'^'' — A/rT4DAoc|7wf 
and D = {ip+\(p-) is the dressing factor . ^^'^^ The right- 
hand-side of the above equation represents the energy 
gain due to the deformation of the phonon state (the 
squeezed effect). This energy gain depends on the spec- 
tral function as well as the phonon density of states 
which is assumed to be p{iu) = u}^~^ /lj^, where n is 
the dimension degree of the bath in the long-wavelength 
approximation.^^ In the Ohmic case (s = 1), it can be 
found that 



ac 



1, for Aq/ud < 1. 



(18) 
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FIG. 2: (Color online) Variation of both the variational 
ground state energy Ei{K") and the displaced oscillator 
ground state energy Eo{K\ = 0) as a function of ct in the 
sub-Ohmic case with a = 0.5 and Aq/ujd = 0.1. Ei(K'i) 
represents the energy maximum lies between. For a > ac, 
Ei{Ki) > Eq{Ki = 0) demonstrates the failure of the varia- 
tional scheme. 



FIG. 3: (Color online) Variation of both the energy difference 
{Eo — El) (with scale on the left) and the non-zero solution 
K'l (with scale on the right. Note that this K'l is the non- 
zero solution in the Ohmic case, not one of the solutions in 
the sub-Ohmic case mentioned above) of self-consistent equa- 
tion with the dissipation strength a in the Ohmic case with 
Aq/ujd ~ 0.1. As a increases, both (Eq — Ei) and K'l decrease 
continuously to zero as q — > 1. 



Once again this result coincides with the crossover 
point predicted by the corresponding variational 
calculation ji^i^ The calculation in the hybrid variational 
scheme is rather complicated and we find it is hard to find 
the crossover point analytically in this case. 

The above analysis provides a different explanation for 
the crossover predicted by the variational calculations. 
When a < ac-, we have Ey < Eq and the variational 
ground state is stable. However, for a > the varia- 
tional scheme fails to give a lower energy and the vari- 
ational ground state becomes unstable, which results in 
a transition from the variational ground state to the dis- 
placed oscillator ground state at a = ac- Since the vari- 
ational ground state has a non-zero tunneling splitting 
while the displaced-oscillator ground state has zero tun- 
neling splitting, such a transition is misidentified as the 
crossover from the delocalized to localized phases, lead- 
ing to an artifact of phonon-induced localization. It is 
important to note that the transition point from the vari- 
ational ground states to the displaced-oscillator ground 
state is determined by Eq. (14). In the Ohmic case with 
Ao/ti'D <C 1, the transition is continuous, while in the 
sub-Ohmic case or Ohmic case with a large Aq/lod, the 
transition is discontinuous and in this case, the transition 
point is given by Eq. (14) according to Landau theoryf^^ 
Based on the above result, it is easy to understand why 
different variational schemes predicted different crossover 
points. 



B. Fail point of the variational scheme at T 7^ 

In the case of T 7^ 0, there is only one kind of varia- 
tional scheme based on the DOS. The application of other 
variational schemes mentioned above is in question since 
the variational scheme based on the DSS was shown to 
fail in thermodynamical limit 1^ Variational calculation 
at T 7^ is different from the case of T = Firstly, 
we perform a displaced-oscillator transformation to the 
Hamiltonian given in Eq. (1), 

H' = U-^HU = Ho + V, (19) 

where 

C/ = expK^/fc(6fe-6i), (20) 

k 

and fk is the variational parameter. It can be found that 
Ho = yfTx + J2^blbk + ^)uJk + (y'lY^iuokfk - 2fkCk), 

k k 

(21) 

V = V+a+ + V-a^ +VQa,, (22) 

where 

V+ = VI = ^exp{-2 5^ A(6fc - bi)} - |, (23) 

k 

^0 = + bk)ick — fk), and A is the tunneling split- 

ting in the presence of dissipation at T 7^ 0. Up to the 
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second order approximation, the up-bound of the free 
energy can be found by the BogoHubov-Feynman theo- 
rem and the result is (setting the Boltzmann constant 

Ab = -Thr[2cosh(/3A/2)] + Y,{ukfl - 2/fcCfe), (24) 



where (3 = 1/T and the variational condition in the 
present case reads 



a = 0.01 , s = 0.45, = OWZa^ 



SAe/Sfk = 0, 



which leads to 

fk = 



Ck 



ujk + A coth(/3wfc/2) tanh(/3A/2) ' 



(25) 



(26) 



and a self-consistent equation of the dressing factor D = 
A/Aq can be found 



D — exp 



-a 



J{Lu)coth{(3u/2)duj 



lo [uj + DAo coth(/3tj/2) tanh(/3A/2)]2 

(27) 

substituting the above result back to Eq. (24) and the 
up-bound of the free energy for the variational thermo- 
dynamical ground state can be found. In analog to the 
case at T = 0, we have a displaced-oscillator thermody- 
namical ground state and its up-bound of free energy can 
be found by setting fk — Ck/wk, the result is 



i^o = -rin2-^ Tdc.:^ 

2 Jo uj 



(28) 



It is clear that the displaced-oscillator ground state corre- 
sponds to the trivial solution of the self-consistent equa- 
tion (27), e.g., Ab{D = 0) = Fq. By the way, as in the 
case of T = 0, we should have < Fq for a successful 
variational scheme and the fail point equation is 



Fn = Af 



(29) 



Numerical analysis shows that, for a fixed temperature, 
the variational scheme fails when a increases to some crit- 
ical point Uc, while for a fixed a, the variational scheme 
fails when T increases to some critical point T^. Typical 
results are shown in Fig. 4. As a matter of fact, the sit- 
uation is more or less that same as in the case of T = 0. 
The variational thermodynamical ground state becomes 
unstable when the fail point is reached and a transition 
from the variational thermodynamical ground state to 
the displaced-oscillator thermodynamical ground state 
happens, which is misidcntified as the crossover from the 
delocalized to localized phases. It is worth noting that 
such a transition happens for < s < 2 in the case of 
T ^ 0.^^ Based on the above analysis, the transition for 
1 < s < 2 (i.e., in the super-Ohmic case) is not a true 
delocalized-localized phase transition, hence it is not in 
confliction with the known result. 
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FIG. 4: (Color online) Comparison of the free energies ob- 
tained from the self-consistent equation in the case of s = 0.45 
and Ao/ti;D = 0.002. The trivial solution _D = represents 
the displaced-oscillator thermodynamical ground state. Two 
non-zero solutions D2 and D\ represents the variational ther- 
modynamical ground state and a maximum lies between, (a) 
Illusion for the failure of the variational scheme for a fixed 
a — 0.01. At the beginning of the low temperature region, 
Ab{D2) < Fo — Ab{D = 0) and the variational thermo- 
dynamical ground state is stable, when T increases to some 
critical point Tc we have Ab{D2) = Ab{D = 0). For T > Tc 
we have Ab{D2) > Ab{D = 0), showing the failure of vari- 
ational scheme at T = Tc. (b) Illusion for the failure of the 
variational scheme for a fixed T/ujd ~ 10~^ in the case of 
Aq/ujd ~ 0.002. The situation is the same as in (a), varia- 
tional scheme fails when q > Qc. 



III. FURTHER ANALYSIS ON FAIL POINTS OF 
THE VARIATIONAL SCHEMES 

In this section, we shall try to answer two questions: 
why the variational schemes fail and why we have differ- 
ent transition behaviors (continuous and discontinuous) 
when the fail point is reached? To this goal, we turn 
to study the oversimplified spin-boson model, i.e., a two- 
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level system coupled to just one phonon mode, the Hamil- 
tonian in the single phonon mode case is given hy^^ 

Hi = -^Aoa^ +u;b^b + Xa:,{P + b). (30) 

The reason why we study this simple model is as follows. 
Firstly, the fail point analysis presented above can be 
used to this model to see whether fail point exists in 
the single mode case. Secondly, the ground state of this 
model can be known analytically in an approximate way 
or more exactly by the numerical diagoanlization^i^i The 
comparison between the variational ground state with 
fail point and the true ground state can help to give the 
answer. 



A. Fail point in the single mode case 

In the case of T = 0, the fail point analysis for multi- 
mode case can be directly generalized to the present case. 
For example, the fail point equation for this model can be 
directly deduced from the corresponding equations given 
before. For clarity, we shall use the variational scheme 
based on the DOS as an example to show the details. 
From Eq. (15), one can find that the fail point equation 
in the present case is given by 



1 



-(fci - fcn)Ao = (Xc/ujy 



(1 + hAo/ioy 



where ki is given by the self-consistent equation 



ki = exp < — 



2A2 



[lo + fciAo)2 



(31) 



(32) 



and kg = e 



. The above equation, in general, have 



three solutions k()*,k[, and /c" where fep* < fcj^ < k'{. Nu- 
merical analysis shows there exists a frequency threshold 
ujc, variational schemes fail only for lu < cuc a-t X = Xc, 
which is w-dependent. It turns out that, for lu > lOc, the 
variational ground state is always stable, i.e., E(ki) < 
E{ko), and the tunneling splitting /ciAq decreases con- 
tinuously with increasing A. However, for lu < iVc, the 
tunneling splitting shows a discontinuous jump as A ap- 
proaching Ac, at where the variational ground state be- 
comes unstable. Figure 5 shows the variation of both the 
energy and tunneling splitting around the fail point. In 
fact, such a discontinuous jump was firstly reported in 
Ref. [3 and stated it as an artifact of the variational cal- 
culation. Our analysis shows clearly this artifact is due 
to the fail point of the variational scheme. Nevertheless, 
the discontinuous jump in the single mode case has some 
important difference from the discontinuous transition in 
the multi-mode case. First of all, the discontinuous jump 
does not exactly lead to a crossover from the delocalized 
to localized phases since both phases have non-zero tun- 
neling splitting. The most important difference is that 
the transition is not from the variational ground state 



to the displaced-oscillator ground state, but to a state 
closed to it. This is because fcg is not a solution of the 
self-consistent equation (32), a situation that is qualita- 
tively different from the multi-mode case where Kq = 
is the trivial solution of self-consistent equation (12) in 
both the Ohmic and sub-Ohmic cases. Since the tran- 
sition determined from the self-consistent equation can 
only happens between the solutions of this equation, the 
transition cannot be from the variational ground state to 
the displaced-oscillator ground state directly, neverthe- 
less the difference becomes very small when oj/Aq < 1/5. 



/ CO = 3.0 




Ap / CO = 3.0 




FIG. 5: (Color online) Illusion of the fail point of the vari- 
ational scheme based on the DOS in the single-mode case 
with Ao /lu = 3.0. In the present case, the self-consistent 
equation has three solutions, k" , kg, and k'l, which represent 
respectively the variational ground state, a state closed to 
displaced-oscillator ground state, and the energy maximum 
lies between, (a) Variation of the two ground state energies 
with A. The variational ground state energy becomes larger 
than the displaced-oscillator ground state when A > Ac, show- 
ing the failure of the variational scheme, (b) Solutions of the 
self-consistent equation as a function of A. When the fail 
point is reached, the variational ground state become unsta- 
ble and tunneling splitting will have discontinuous jump, i.e., 
from fc" Ao to kgAo- 

Fail points of other variational schemes can be ana- 
lyzed in the same way. The situation is almost the same 
for all three variational schemes except the values of cuc 
and the corresponding Ac. Figure 6 shows the boundary 
lines of the fail points as a function of Aq /lu and X/lu for 
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all the three variational schemes. As one can see from 
Fig. 6, all variational schemes fail in the low frequency 
and strong coupling regions. 
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FIG. 6: (Color online) Fail point boundary lines of the three 
variational schemes in the single-mode case. The starting 
points of the curves are (Ao/a;c — 2.54; Ac/cj — 1.37), 
(Ao/cJc ^ 4.67; Ac/tJ ~ 1.70), and (Ao/a;^ ~ 2.46; Kjio ~ 
1.61) respectively for the DOS, DSS, and hybrid case. 

It should be noted that the above analysis cannot be 
generalized directly to the case of T ^ 0. This is because 
the system with a Hamiltonian given in Eq. (30) is not 
a true thermodynamical system. To perform the varia- 
tional calculation as in the multi-mode case, one needs 
to include the environment as a "thermo-bath" , i.e., the 
effect of the environment is just to keep the single- mode 
system in thermodynamical equilibrium. Under this 
sense, we can do a variational calculation based on the 
DOS and the up-bound of the free energy in the present 
case is given by 

U = -Tln[2cosh(/5fc2Ao/2)] -tj-i(C2 - 2CA), (33) 

where 

^ 1 + k2AoUJ-^ coth{/3uj/2) tanh(/3fc2Ao/2) ' ^^"^^ 

and k2 is determined by the self-consistent equation 

_ 2A2coth(/3w/2) 
'^^^[uj + fczAo coth(/3tj/2) tanh(/3fc2Ao/2)]2 ' 

(35) 

Correspondingly, the up-bound of the free energy of the 
displaced-oscillator ground state is 

/o - -T ln[2 cosh(/3A:3 Ao/2)] - X^lo, (36) 

where — exp{— 2(A/cj)2 coth(/3a'/2)}. Using the above 
result, one can make analysis on the stability of varia- 
tional thermodynamical ground state as before. As ex- 
pected, it turns out that the situation is almost the same 



as in the case of T = 0. There is a frequency threshold lOc- 
For uj > ujc, we have fa < fo, i.e., variational thermody- 
namical ground state is always stable for all the coupling 
range and fc2Ao decreases continuously with increasing 
A; however, for lo < Uc, variational thermodynamical 
ground state becomes unstable as the coupling reaches 
some critical value Ac which is again w-dependent. Tun- 
neling splitting shows a discontinuous jump when the fail 
point is reached. In the case of T 7^ 0, the fail point can 
be reached in another way; that is, for a given coupling A, 
one can find that the variational scheme fails when tem- 
perature rises to some critical value Tc, at where tun- 
neling splitting also shows a discontinuous jump. The 
boundary line of the fail point for the variational scheme 
based on the DOS in the case of T ^ is shown in Fig. 7. 
The boundary lowers down as the temperature increases, 
which implies that the variational scheme fails at a lower 
Ac as T increases. 
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FIG. 7: (Color online) Boundary lines of fail point for vari- 
ational scheme based on the DOS in the single-mode case at 
T 7^ 0. As T raises, the boundary lines lowers down, showing 
that the variational scheme fails at a lower Ac as T increases. 

Why variational schemes fail in the low frequency re- 
gion? Physically, in the high frequency region with 
[jj/Aq ^ 1, the first term in Eq. (1) or (30) can be treated 
as a small perturbation, hence the DOS is a good trial 
ground state for the variational calculation; however, in 
the low frequency region with Ao/ti^ ^ 1, the DOS is not 
appropriate to serve as the trial ground state since the 
term iAoCTa; makes the main contribution in the single 
mode case. As one can see from the analysis presented in 
Ref. 4, in the low frequency region, the true ground state 
of the phonon bath in the single mode case is a combi- 
nation of the two DOS with opposite displacement, i.e., 

= ^=^[e-^'i(^^"-^) +ae^'^(^'-^)]|0), (37) 
V 1 + a 

where g\^2 and a are functions of A/Aq and w/Aq. Nu- 
merical analysis shows that \'^ph) gives a very accurate 
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description for the phonon ground state in the low fre- 
quency region with strong couphngs. Obviously, \(l)ph) 
can be described by neither a DOS-like nor a DSS-like 
state, showing that both the DOS and DSS can not serve 
as the trial ground states in the low frequency region 
with strong couplings. As shown in Ref. 4, a variational 
scheme by taking a trial ground state as \(t>ph) shows no 
discontinuous jump. We have made numerical analysis on 
the energy of the variational ground state based on \(f>ph), 
which is found to be always stable against the displaced- 
oscillator ground state in all the coupling regime. This 
result clearly shows that the failure of the variational 
scheme is due to the invalid trial ground state in the low 
frequency regions. The above analysis provides the an- 
swer to why the variational schemes based on both the 
DOS and DSS always fail in low frequency region with 
strong couplings. 

It should be noted that the low frequency modes play 
important roles in the localized transition. As one can 
see from the self-consistent equations (12) or (27), the 
integral can be divided into two parts like, /g ° — s- 
lo'^ -l- J^" , without the contribution of the low frequency 
modes lie in (0, Wc), there is no localization for any finite 
a. In other words, it is the low frequency modes lead 
to localization. The failure of the variational schemes 
in the low frequency regions implies that the variational 
calculations fail to take account of the contribution of 
the low frequency modes to the localization, making the 
predicted localization in question. 



B. Continuous or discontinuous crossover 

To the first order approximation (i.e., by omitting the 
cooperative effect between different phonon modes), the 
result of multi-mode case can be considered as a com- 
bination of all the phonon modes lie in {0,lud) with a 
frequency-dependent weight. Generally speaking, the 
parameter s in the spectral function J{uj) determines 
the weight, while for counting the contribution of the 
phonon modes that have fail point; the parameter Aq/lojj 
is therefore important and the boundary lOc is in scale 
of Aq as shown in Fig. 6. Let us take the variational 
scheme based on the DOS as an example. Giving that 
/S.q/loo = 0.1, one can see from Fig. 6 that Wc/Aq ~ 1/3, 
which implies that the frequency modes with frequency 
< a; < O.OSwd have fail points, i.e., about 3 percent 
phonon modes have fail points. However, if one takes 
Ao/w£) = 1 as shown in Fig. 1, then there are about 
30 percent phonon modes have fail points. This result 
shows that the population of the phonon modes having 
fail points is directly controlled by the value of 
Now we know that only the phonon modes have fail 
points will contribute a discontinuous jump of the tun- 
neling splitting as the coupling reached the critical point. 
In the Ohmic case with Aq/wij <C 1, the discontinuous 
contribution of the low frequency modes at the fail point 
is overwhelmed by the continuous contribution from the 



high frequency modes, leading to a continuous transition; 
however, as Aq/wu increases to some larger value (say, 
0.5 in the DOS case), the increasing population of phonon 
modes having fail points leads to an observable discon- 
tinuous contribution and the transition behavior turns 
to discontinuous. In the sub-Ohmic case, the weight 
of the low frequency modes increases so that the dis- 
continuous contribution from phonon modes having fail 
points becomes observable even when Aq/wd 1, hence 
the transition is always discontinuous. However, for the 
same reason, one can easily predict that the discontinu- 
ous transitions will be unapparent when Aq/w^) becomes 
too small. It should be noted that the discontinuous tran- 
sition in the sub-Ohmic case is just due to the fail point 
of the variational schemes, not the true crossover behav- 
ior of the model since it is not expected to be seen in the 
NRG or DMRG calculation.^^i25 

The above analysis can be generalized to the case of 
T 7^ 0. It has been shown that the transition behav- 
ior predicted by the variational calculations at T 
is discontinuous for < s < 2i2^ The key point is to 
understand why the discontinuous transition behavior is 
extended to the super-Ohmic case. By comparing the 
two self-consistent equations (12) for T = and (27) for 
T 7^ 0, at T 7^ 0, the spectral function is effectively mod- 
ified as J(w) coth(/3a'/2), i.e., an extra factor coth(/3cL'/2) 
appears in numerator of the integrand. It is obvious that, 
coth(/3cj/2) oc 1/uj as cj ^ 0, which implies that, in the 
low frequency regions, J{uj) coth(/3a;/2) cx In other 

words, the extra factor coth(/3ti;/2) seriously increases the 
weight of the low frequency modes having fail points, re- 
sulting in a discontinuous transition behavior even in the 
super-Ohmic case. 

Mathematical analysis on the self-consistent equations 
derived from the variational condition at both T = 
and T shows there exists a universal Sc for all the 
variational schemes, the predicted transition is always 
discontinuous for < s < Sc, a continuous transition can 
only be found when s — Sc with Aq /lud <C 1. It was 
found that = 1 for T = and = 2 for T 7^ 0.^^ ^he 
above analysis clearly shows that the predicted transition 
behavior depends on both the weight and the population 
of the low frequency modes. 



IV. CONCLUSION AND DISCUSSION 

In conclusion, we have shown that all the variational 
schemes used in the spin-boson model have fail points. As 
the fail point is reached, the variational schemes fail to 
give a lower energy (or lower free energy at T 7^ 0) than 
the displaced-oscillator ground state, leading to a transi- 
tion from the variational ground state to the displaced- 
oscillator ground state, which is always misidentified as 
the crossover from the delocalized to localized phases in 
a variational calculation. Our analysis show that the ex- 
isted variational schemes are not suitable for studying the 
crossover problem in the spin-boson model. The present 
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analysis can help to understand the controversy between 
different variational schemes and discrepancy with other 
calculations, like the NRG and DMRG. The analysis on 
the fail point in the single mode case can help to answer 
why the variational schemes have fail points and how dif- 
ferent transition behaviors were found in different condi- 
tions. However, the situation is still far from totally sat- 
isfactory. One intriguing fact is that the "fake" crossover 
point predicted by the variational calculation based on 
the DOS is found to be in good agreement with the re- 
sults obtained by other methods including the NRG cal- 
culation for both the Ohmic and sub-Ohmic cases i^'^^i^^ 
In the present stage, it is not yet clear if this is just a 
coincidence. One possible explanation is that since the 
population of the low frequency modes having fail points 
is very small in the case of Ao/wd <C 1, the contribution 
of them is small comparing with the rest of the phonon 
bath. Our DMRG calculation also assures that the vari- 
ational ground state based on the DOS is a good approx- 
imation to the true ground state before the transition 



happens."^ Nevertheless, as we have mentioned before, 
the phonon-induced localization is mainly due to the low 
frequency modes, combined with the discontinuous tran- 
sition behavior in the sub-Ohmic case, the answer is not 
very convincing. Perhaps one can find that answer by 
developing a variational scheme without fail point and 
this work is now under process. The present result can 
also shed some lights on the failure of the variational 
calculation on the self-trapped transition in ploaron and 
exciton systems. In fact, the present analysis is read- 
ily applicable to the discontinuous transition predicted 
by the variational calculation gave in Ref. 4, but further 
analyses are needed for other cases i^^^i 
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